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ABSTRACT 



Aims. We derive rotation curves from optical emission lines of 182 disk galaxies (96 in the cluster and 86 in the field) in the region of 
Abell 901/902 located at z ~ 0.165. We continue the kinematic analysis presented in a previous paper. Here, we focus on the analysis 
of B-band and stellar-mass Tully-Fisher relations. We examine possible environmental dependencies and differences between normal 
spirals and "dusty red" galaxies, i.e. disk galaxies that have red colors due to relatively low star formation rates. 
Methods. Assuming a tilted ring model we simulate the spectroscopy of a given galaxy and reconstruct the maximum rotation ve- 
locity from the best-fitting simulated rotation curve. We fit regression lines adopting a maximum likelihood method based on Bayes' 
theorem to build B-band and stellar-mass Tully-Fisher relations. 

Results. We find no significant differences between the best-fit TF slope of cluster and field galaxies. At fixed slope, the field pop- 
ulation with high-quality rotation curves (57 objects) is brighter by AM B = -0M2 ± 0M5 than the cluster population (55 objects). 
We show that this slight difference is at least in part an environmental effect. The scatter of the cluster TFR increases for galaxies 
closer to the core region, also indicating an environmental effect. Interestingly, dusty red galaxies become fainter towards the core at 
given rotation velocity (i.e. total mass). This indicates that the star formation in these galaxies is in the process of being quenched. 
The luminosities of normal spiral galaxies are slightly higher at fixed rotation velocity for smaller cluster-centric radii. Probably these 
galaxies are gas-rich (compared to the dusty red population) and the onset of ram-pressure stripping increases their star-formation 
rates. 

Galaxies of smooth morphology show higher rms values in the fitting of their rotation curves. This is particularly the case for dusty 
red galaxies. A cluster-specific interaction process like ram-pressure stripping is the best explanation, since it only affects the gaseous 
disk but not the stellar morphology. 

Conclusions. The results from the TF analysis are consistent with and complement our previous findings. Dusty red galaxies might 
be an intermediate stage in the transformation of infalling field spiral galaxies into cluster SOs, and this might explain the well-known 
increase of the SO fraction in galaxy clusters with cosmic time. 

Key words, galaxies: clusters: general - galaxies: clusters: individual (A901, A902) - galaxies: evolution - galaxies: kinematics 



1. Introduction 

The empirically found, tight correlation between the rotation ve- 
locity V mt of disk galaxi es and their luminos ity L(L cc V" ot ) was 
originally introduced by Tully & Fisher] (ll977l) as a new method 
of determining distances to galaxies. Since the maximum rota- 
tion velocity is proportional to the total mass of a galaxy and the 
luminosity correlates with stellar mass, the Tully-Fisher relation 
(TFR) provides a link between the dark and visible matter com- 
ponent of a galaxy. This link is assumed to be a consequence of 
the formation scenario of di sk galaxies withi n the framework of 
a cold dark matter universe (iBlumenthal et a l. 1986). Initially in 
equilibrium with the dark matter halo, rotating disks form from 
collapsing gas clouds. In the process, angular momentum is con- 
served. In contrast to dissipationless dark matter, baryonic par- 
ticles can loose part of their energy via radiation and hence are 
able to form more compact structures (e.g. Dutton et al. 2007). 



Analysing the Tully-Fisher relation in different luminosity bands 
(e.g. Kannarjpan et al. 2002; Pizagno et al. 2007; Verheiien 
2001) or plotting the rotati on velocity d irectly against the 
stellar mass (smTFR; e.g . iKassin et ail l2007t iMiller et all 
2011; Pizagn o et al.l l2005|) or baryonic mass (b TFR; e.g. 
McGaugh *et all I2000t1puech et all [20lok IVerheijerJ l200ll) can 
give valuable insight into the formation history of spiral galaxies. 
Comparing the best-fit parameters of the TFR (slope, intercept 
and scatter) at different redshifts to predictions of numerical sim- 
ulations helps to constrain the physica l processes important for 



the evolution of spiral galaxies (e.g. Bamford et al. 2006, 2005 
Bohm & Ziegle^|2003 'Bohm et al. 2004; IKassin etalj 12 007; 
Mill er et alJl201 ll 120121: lPuech et al.ll2008TlVogt et al.lll996l). 



* Based on observations with the European Southern Observatory 
Very Large Telescope (ESO-VLT), observing run ID 384.A-0813. 
** e-mail: benjamin.boesch@uibk.ac.at 



The well-known morphology-density relation (Dressier 1980) 
has been confirmed in many subsequent studies highlighting a 
picture of high-density environments altering the properties of 
galaxies. Besides a prevalence of more early-type morphology, 
galax ies residing in cluster s on average have e.g. redder colours 
(e.g. Blanton et al. 2005), lower star-formation rates (e.g. 
[Lewis et alJ 120021 : IVerdugo et al.l l2008) and a reduced gas con- 
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tent (e.g. iGiovanelli & Havnel 19851: ISolanes et alJl200lh com- 
pared to field galaxies. Consequently the Tully-Fisher parame- 
ters might vary between different environments and hence are 
the topic of several stu dies (e .g. Bamford et al. 2005; Jager et al. 
2004; Kutdemir et all 120101 12 008; Milvan g- Jensen et al J 12 003; 
Moran et al. 2007; Nakamura et al. 2006; Ziegler et al. 2003). 
Some of them show conflicting results and conclusions. The 
main reason for this are differences in the data analysis includ- 
ing the observation of the target galaxies, the derivation of lu- 
minosities and rotation velocities, or the methods used to fit the 
TFR. Several definitions of rotation velocities exist and many 
groups use slightly different methods to account for observa- 
tional effects and to infer the intrinsic rotation curve of a galaxy. 
Different sample sizes, selection criteria and galaxy populations 
make a comparison of results even more complicated. 
However, several general trends could be established. The slope 
of the TFR increases towards redder photometric passbands (e.g. 
iKannappan et alj 120021; iPizagno et all 120071 ; IVerheiierj 1200 lh . 
Redder bands are less sensitive to recent star-formation events 
and dust extinction and better trace the underlying old stel- 
lar population, which dominates the stellar mass. On the one 
hand, the B-band TFR shows a brightening of the intercept 
up to AM B ~ 1 mag out to a redshift of z ~ 1 (e.g 
Boh m et all 12004 "Fernande z Lorenzo et alj 1201 Ot [M iller et al 



201l UVogtetaijT l996: Wei ner et al.ll2006UZiegler et al.l l2002 



This is most likely due to an enhancement of the sta r-form ation 
rate and a younger stellar population (Hopkins 2004]). On 
the other hand, there is little or no evolut ion in the stellar- 
mass or baryonic Tully-Fisher relati on (e.g. iFlores et all 12006; 
Miller et al. 201 1; Puech et al. 2008). In agreement with simu- 
lations (e.g. [Steinmetz & Navarro 1999), most s tudies measure 
a larger scatter for TFRs at higher redshifts (e.g . Bamford et al.l 
l2006llBohm & Zieglerjl2007tlKassin et alll200% . 
Comparisons between the TFRs of cluster and field galaxies 
show no cl ear results. Wh ile Jaffe et all (1201 lh , [Nakamura et al.l 
(2006) and Zi egler et alj d2003l) fin d no significan t differ ences 
between these two environ ments, iBamford et al.1 ([2005) and 
Milvang- Jensen et al. (2003) deduce higher B-band luminosities 
in clusters compared to the field. Moran et al . ( 2007 ) find a larger 
scatter for cluster galaxies. Again, sample selection and analy- 
sis techniques might influence the results. But also the choice of 
a particular cluster and its ICM properties and dynamical state 
might play a role. 

Using spatially resolved spectra obtained with VIMOS (Visible 
MultiObject Spectrograph) on the VLT (Very Large Telescope) 
in combination with already existing data from the STAGES 
(Space Telescope A 90 1/902 Galaxy Evolution Survey) project 
(Grav et al. 2009a), we continue the kinematic analysis of 
Bosch et al. (2013, hereafter Paper I). We exploit rotation curves 
of disk galaxies in and around the multiple cluster system Abell 
901/902 with a focus on the analysis of the Tully-Fisher relation. 
The paper is organised as follows. Sect. 2 gives a short introduc- 
tion into the STAGES project, outlines the sample selection and 
observations and reviews results of Paper I. Sect. 3 describes the 
derivation of the parameters needed for the Tully-Fisher analysis 
and presents adopted methods to build the TFR. Sect. 4 contains 
the TF analysis, which provides the basis for the residual analy- 
sis of Sect. 5. Finally, Sect. 6 summarises the main points. 
Throughout this paper, we assume a cosmology with Q m = 0.3, 
Qa = 0.7 and Ho - 70 kms^'Mpc -1 . Magnitudes quoted in this 
paper are in the AB system. 



2. Data 

2.1. The STAGES data set 

STAGES is a multi-wavelength survey centred on the area 
around the multiple cluster system Abell 901/902. It comprises 
four subclusters: A901a (at redshift 0.164), A901b (z=0.163), 
A902 (z=0.167) and the SW group (z =0.169). They are in 
a non-virialised state dBosch et al.ll2013l) most likely coalesc- 
ing to one big system in the future. The 0°5 x 0°5 (~ 5 x 
5Mpc 2 ) region of Abell 901/902 has been the subject of Hubble 
Space Telescope/ Advanced Camera for Surveys (HST/ACS) ob- 
servations producing an 80-tile mosaic of V-band (F606W) 
images. We have at hand 17-band COMBO-17 (Classifying 
Objects by Medium-Band Observations) ground-based imag- 
ing d Wolf et al J l2003h. Spitzer 24 /urn, and XMM -Newton 
X-ray data (IGi lmour et al. 2007). Gravitational len sing maps 



( Hevmans et al.ll2008l) . spectral energy di stributions dWolf et al.l 
2003), stellar masses dBorch et al.l 120060 . star-formation rates 
(Belletal. 2007), mean stellar ages and morphological classi- 
fications provi de an e xcellent starting point for further analy- 
sis. [Gra^eTa]] (|2009a|) give an overview of the STAGES project 
and pr esent the publicly available master catalogue dGrav et al.l 
2009B). 



2.2. Observation, sample selection & data reduction 

Information on the spectroscopic observations of A90 1/902 with 
VLT/VIMOS as well as details on the data reduction process are 
presented in Paper I. Thus, we will only give a short review of 
the most important steps here. 

Multi-object spectroscopy (MOS) with VIMOS was carried out 
between February 8 and March 10, 2010 (ESO-ID 384.A-0813, 
P.I. A. Bohm). Four VIMOS pointings (~ 1' overlap) sample 
the whole A90 1/902 field (two exposures of 1800 seconds 
each). We opted for the high resolution grism HR-blue with a 
slit width of 1'.'2, which covered a ~2050A spectral range in 
the wavelength interval [3700A, 6740A]. This configuration 
yielded a spectral resolution of R ~ 2000 and an average 
dispersion of 0.51 A/pixel at an image scale of 0V205 /pixel. 
The robustness of our analysis against varying seeing conditions 
during spectroscopy are discussed in Sect. 13.31 
Our primary targets for spectroscopy are disk galaxies showing 
emission lines. Hence, we selected galaxies with stellar masses 
M„ > 1O 9 M , absolute magnitudes Mb < -18, a s tar-forming 
SED (selection criterion in the public catalogue iGrav et akl 
(2009b): sed_type > 2) and a visually confirmed disk compo- 
nent on the ACS images with an inclination angle i > 30°. We 
require photometric redshifts in the interval 0.1 < z p hot < 0.26 
resulting in a matched cluster and field sample. This selection 
resulted in ~320 cluster and ~160 field candidates. Using the 
ESO software VMMPS (VIMOS Mask Preparation Software), 
we manually placed MOS slits along the apparent major axis of 
the galaxies (allowing slit tilt angles up to +45°). In the end, we 
obtained spectra of 215 different galaxies. 
One science driver of our project is to investigate the population 
of dusty red galaxies (Wolf et al. 2005). In total, 44 of the 215 
spectra belong to galaxies with this SED ty pe (selection crite- 
rion in the public catalogue IGrav et all (l2009bl) : sed_type = 2). 
Dusty red galaxies are spread across the red sequence in the 
U-V colour-magnitude diagram. But, unlike early-type red 
galaxies, part of their colour is due to intrinsic dust extinction 
(E(B - V) > 0.1). Hence, they co uld be interpreted as the low 
specific-SFR tail of the blue cloud dWolf et al.ll2005h . 
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Fig. 1. Redshift histogram of the 200 galaxies of our sample. The 
cluster sample (106 objects; blue, "backslash" hatching) peaks 
at the cluster redshift z ~ 0.165. The field sample (94 objects; 
green, "slash" hatching) shows a bump at z ~ 0.26, which might 
be a volume effect. 



We performed the spectroscopic data reduction mainly using 
the ESO-VIMOS pipeline (version 2.2.3). The main reduction 
steps were bias subtraction, correction of distortions resulting 
from the focal reducer using flat field exposures (a potential 
normalisation was skipped due to reflections present in the 
flat-field images) and wavelength calibration. Additionally, 
we performed bad pixel and cosmic ray cleaning. Wavelength 
calibration for spectra with slits in the upper third of the CCD 
or with large tilt angles was improved by adding Argon lines to 
the catalogue of arc lamp lines. For more details see Paper I, 



2.3. Results from Paper I 

In thi s subsection we will review results from iBosch et all 
(120131) . which will be important for the current analysis and help 
to interpret our data. 

We were able to measure spectroscopic redshifts for 200 dif- 
ferent galaxies, either from emission lines (188) or from strong 
absorption lines (12). Using galaxies obser ved more than once 
and fo llowing a prescription outlined in iMilvanp- Jensen et al.l 
(2008), we estimated the typical redshift error to be 6 Z /(1 + z) ~ 
0.00016 or 47km/s in rest-frame. 

We used the biweight estimator of scale and location 
dBeers et al.ll 19901) to estimate the velocity dispersion and mean 
redshift of each of the four subclusters individually. On the basis 
of these properties, we determined cluster membership using a 
3-cr-clipping. We identified 106 cluster members. Fig. Q]shows 
the redshift histogram of the whole sample, separated into clus- 
ter and field galaxies. 

We found that the fraction of cluster galaxies with distorted 
rotation curves is 75% higher than for field galaxies. Most of 
these differences can be attributed to morphologically undis- 
torted galaxies with a relatively smooth stellar disk. Dusty red 
galaxies account for a large part of this population. This suggests 
a cluster-specific mechanism, which mainly affects the gas disk 
of the cluster spirals, with ram-pressure stripping as the best can- 
didate. Investigations of the gas concentration implied that most 
spiral galaxies traverse the cluster for the first time. Especially 



for dusty red galaxies, ram-pressure could be a dominant mech- 
anism in A901/902, even in the cluster outskirts. 



3. Methods 

Most galaxy properties such as stellar masses, apparent mag- 
nitudes and inclination angl es are already provid ed by the 
STAG ES project (Gray et al. 2009bj ). Assuming a Kroupa et al.l 
(Il993l) stellar initial mass function. iBorch et aD (2006) derived 
stellar masses from an S ED-fitting of templates from PEGASE 
(Fioc & Rocca-Volmerange 1997) population synthesis models 
to the COMBO- 17 fluxes. They noted that alternatively using 
Kroupa (2001) or Chabrier (2003) initial mass functions would 
result in the same stellar masses to within ~ 10%. 
We only need to determine rest-frame absolute magnitudes and 
maximum rotation velocities for each object, and then derive the 
Tully-Fisher relations for several samples. We adopt V2.2, the ve- 
locity at 2.2 disk scalelengths r,/, as our kinematic measure, be- 
cause it produces the tightest TFRs and is more robu st against 
observational errors than other velocity definitions (ICourteaul 
1997). In contrast, the asymptotic velocity V max , e.g., is easily 
overestimated in galaxies with insufficiently extended rotation 
curves. We note, that the chosen point at 2.2r</ is the radius where 
the rotation vel ocity of an ide alised self-gravitating disk is ex- 
pected to peak (iFreemanll 19701) . 



3.1. Rest-frame magnitudes 

Apparent V-band magnitudes MagJjest from the HST/ACS 
F60 6W) images we re measured using SExtractor 
Berli n & Arnoutsl I1996I) . To compute absolute magnitudes 
in the B-band, we apply a k-correction ks(V, T,z), the distance 
modulus D(z, Ho, O m , Q./0 and a correction for Galactic absorp- 
tion Ay and intrinsic extinction A' B {a/b, V max ). 
The k-correction is needed to transform the apparent magnitude 
in an observed band to a rest-frame band for a given galaxy. 
This correction depends on the SED type T and the redshift z. In 
our case it is sufficient to use the transformation V — » B, since 
at the redshifts of our sample the V-band matches the rest-frame 
B-band. We determine the correct SED T-type individually by 
comparing our galaxy spectr a to the templates in the Kennicutt 
catalogue (Kennicutt 1992a.b). Besides the stellar continuum 
shape of the spectra, we use e.g. equivalent widths, the relative 
strength of emission lines, and the prominence of absorption 
lines for classification. The SED type T ranges from 1 (for a 
typical spectrum of an Sa galaxy) to 10 (for type Im). For more 
details see lBohm et all (120041) . 

We adopt the correction factors for Galactic extinction in 
the direction of A901/902 (a=9h56minl7s, <5=-10°01'll", 
J2000) of the obse rved V-band from the NASA extragalactic 
database following Sch laflv & Finkheinerl d201ll) . This yields 
the following value for the whole sample: A„ = 0^1 36. 
A critical step is the correction for intrinsic absorption 
A' B (a/b, V max ) due to dust extinction. The intrinsic absorption 
depends on the inclination angle i or equivalently the axis ratio 
a/b of a galaxy and is maximal for an edge-on and mi nimal for a 
face-o n configuration. We follow the prescription of iTullv et al.l 
( 1998) and apply a correction, which is also correlated with the 
maximum rotation velocity V max of a galaxy, s ince extinction 
is stronger in more massive galaxies (Giovanelli et al. 1995). 
This correction is relative to a face-on orientation and does not 
account for a potential residual absorption in a face-on system. 
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All in all, we apply the following intrinsic absorption correction: 

A B (a/b, V max ) = log(a/b) [-4.48 + 2.75 log(V max )] (1) 

To summarise, the absolute magnitude in a galaxy's rest-frame 
B-band is calculated as follows: 

M B = m v - D(z) - k B (V, T, z) - A s v - A B (a/b, V max ) (2) 

We compute the corresponding errors of the absolute magnitudes 
via error propagation. We assume an uncertainty of 10% for the 
axis ratio and an SED classification error of 6T = 2. The lat- 
ter transforms into an average uncertainty in the k-correction of 
cr k = 0TO8. 

3.2. Rotation-curve extraction 

We address the extraction of the rotation curves from the emis- 
sion lines of a spectrum only very briefly here, since this process 
is already described in detail in Paper I. Depending on a galaxy's 
redshift and position on the CCD, up to four emission lines are 
visible in the VIMOS spectra: [OIII] A5007A, [OIII] A4959A, 
H/3A4861A and [Oil] A3! 26/ 31 29. Prior to the emission line 
fitting we add up three (or five, in cases with a very low S/N) 
neighbouring rows for each spatial position of the spectrum to 
enhance the S/N. We then transform red- and blueshifts of the 
emission lines into an observed rotation curve. The kinematic 
centre of a galaxy by definition has no Doppler shift. Its exact po- 
sition significantly influences the shape of the rotation curve and 
must therefore be derived with care. Primarily, we determine the 
kinematic centre of a galaxy to within 07 1 by fitting a Gaussian 
to its (spatial) luminosity profile at the spectral position of the 
emission line, i.e. the luminous centre is assumed to be also the 
kinematic centre. In a second step, we slightly shift the centre 
by up to ±1.5 pixels, minimising the rotation curve asymmetry 
adopting an asymmetry measure defined in Dal e et all (120011) ). 
We extracted rotation curves from at least one emission line of 
182 different galaxies (86 in the field and 96 in the cluster). 

3.3. Rotation-curve modelling 

At face value, an extracted rotation curve (RC) only traces the 
line-of-site rotation velocity of a galaxy. This observed position- 
velocity information does not represent the intrinsic rotation of 
a disk galaxy. Several geometric and observational effects have 
to be taken into account. The inclination of the disk with respect 
to the line-of-sight, the position angle of the galaxy's apparent 
major axis with respect to the direction of dispersion, beam- 
smearing due to the integration along the slit width and seeing 
effects may "disguise" the in trinsic shape of th e rotation curve. 
For a detailed discussion see Boh m et all (|2004). 
Modelling the galaxy as an infinitely thin rotating disk (tilted 
ring model) with an exponential intensity profile and simulating 
the spectroscopic observation process produces a synthetic rota- 
tion curve. Matching the synthetic RC thus obtained with the ac- 
tual observed rotation curve allows to infer from the best-fitting 
simulated RC the true intrinsic rotation-curve parameters . 
As input parameters for the whole simulation process, we use 
the redshift of the galaxy, the FWHM (Full Width at Half 
Maximum) of the seeing disk, the photometric scale length r</, 
the inclination angle i and the position angle 9 of the galaxy. We 
obtain the photometric scale length rj by fitting a 2-component 
model with GALFIT (Pengetal. 2002) to the corresponding 
HST/ACS V-band images. We assume a de Vaucouleurs pro- 
file for the bulge and an exponential profile for the disk. We 



stress that the rotation curve modelling explained below is robust 
against the assumed morphology, and changes in the morpholog- 
ical decomposition (like e.g. a free-n Sersic bulge) do not affect 
our results. Values of the axis ratio ajb (a and b are the apparent 
semi-major and semi-minor axis, respectively) and the position 
angle were previously deter mined using the SExtractor pack- 
age (iBertin & Arnou ts 1996) on the HST/ACS V-band images 
and ar e adopted from the STAGES master catalogue (Gray et al. 
2009b). Assuming an intrinsic disk thickness of q — 0.2, we 
computed the inclination from the axis ratio (e.g. iTullv et all 
1998). 

A robust value for the FWHM of the seeing disk is an impor- 
tant ingredient in correctly simulating the observation process. 
The presence of three reference star spectra in each quadrant 
of a pointing allows to determine the seeing conditions during 
spectroscopy. This is clearly superior to estimates inferred from 
the Differential Image Motion Monitor (DIMM). We determine 
the FWHM of the stellar spectra along the spectral axis for each 
MOS image. We adopt the median of the three reference stars as 
our input parameter for the simulation. Seeing values increase 
towards shorter wavelengths and a standard equation to describe 
this variation is given bv lSarazinl (120031) : 



I A \ 1/5 
FWHM(/l) = FWHM(/tc) I — I . 



(3) 



where we use A c - 5100A as the central wavelength of the HR- 
blue grism. We obtain seeing values between 0760 and 1727 with 
an average of 0793. 

In the following, we describe the simulation process in more de- 
tail. We assume an intrinsic rotation function of the following 
form (ICourteaulll997l) : 



CM = 



Vn 



(r° + r?) 



«\l/o' ' 



(4) 



At the turn-over radius r, the linear rise for small radii (rigid- 
body region) turns over into a flat region with asymptotic ve- 
locity V max , where the dark matter halo dominates the mass dis- 
tribution. The factor a determines the s harpness of the turnover 
and is set to a — 5 (Bohm et al. 2004). Hence, V max and r t are 
the two remaining free parameters in the fitting process. 
As a starting point for the simulation process, we generate a two- 
dimensional (2D) exponential intensity profile on a pixel grid: 



I{x,y) -he 



(5) 



where r - yx 2 + y 2 and h is the intensity at the galactic centre 
at the pixel position (x,y) - (0,0). Then we superimpose (on 
this intensity profile) a 2D velocity field using the above rota- 
tional law (Eq.|4]l. This results in a rotating disk seen face-on. In 
the next step, we obtain the correct geometric projection. To this 
end, we rotate our galaxy, i.e. the 2D intensity and velocity field, 
according to the inclination angle i (rotation around the y-axis) 
and the position angle 8 (rotation around the line-of-sight or z- 
axis). For the velocity field we calculate the line-of-sight com- 
ponent only, since this is the component that is responsible for 
the Doppler shift and hence measured in a spectrum. We show 
an example of such a projected 2D intensity and velocity field in 
Fig. EKleft panel). 

Besides geometric projection effects we also have to consider the 
blurring effect due to seeing. We take this into account by con- 
volving the intensity and velocity field with a Gaussian point 
spread function (PSF) (the FWHM is obtained from the see- 
ing value during spectroscopy). Then, as in the observations, 
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Fig. 2. Simulated 2D intensity and velocity field of a galaxy. The fields are projected according to the inclination and position angle 
(here i - 40°, = 30°). The direction of dispersion is along the y-axis. The exponential intensity profile is depicted in blue colour 
scale (here the scale length is rj = 5.3 pix). Isophotes corresponding to [50, 30, 15]% of the central intensity are superimposed (black 
dotted ellipses). The coloured curves ("spider-diagram") indicate iso-velocity lines ranging from -60%V max (red) to 60%V max (blue) 
in steps of 10%. Right panel: The intensity and velocity field is additionally convolved with a Gaussian Point Spread Function 
(FWHM = 077 = 3.5 pix). The slit (in black) is placed along the apparent major axis of the galaxy. At each spatial position along 
the x-axis, an intensity-weighted average of the velocity values inside the slit is calculated. The shaded magenta area indicates such 
a velocity range for a certain spatial position (with boxcar smoothing of 3 pixels). Here, the kinematic and the luminous centre 
coincide. 



we place a slit along the apparent major axis of the galaxy. 
For the following slit integration the velocity field is weighted 
by the normalised intensity profile, since velocities stemming 
from brighter regions contribute more to the signal in a spec- 
trum. Prior to this intensity weighting of the velocity field, we 
shift the intensity profile according to the mismatch between 
the kinematic and luminous centre, which we determined in the 
course of rotation-curve extraction (see Sect. 13. 2t . Then, for each 
pixel position along the spatial axis (x-axis) the velocity values 
within the slit borders are integrated along the direction of dis- 
persion (y-direction). Note that the integration is in general not 
perpendicular to the slit direction since most MOS slits are tilted. 
Effectively, we calculate an intensity-weighted average of the ve- 
locity values inside the slit at each spatial position. For the inte- 
gration area at each spatial position we account for the size of 
the boxcar filter (3 or 5 pixels) that was used for the extraction 
of the rotation curve. Thus, we obtain a simulated value for the 
rotation velocity at the position of each pixel along the spatial 
axis. This position-velocity information is the end product of our 
simulation and defines a synthetic rotation curve. The right panel 
of Fig. [2] shows a scheme of the integration process of the PSF- 
convolved and intensity-weighted velocity field. Unlike shown 
in the figure, we actually oversample the pixel grid by a typical 
factor of five to improve the precision of the computation. 



3.4. Rotation-curve fitting 



For each galaxy we vary the free parameters (V max , r,) of the 
synthetic RC such that it best fits the observed RC. We do this 



by minimizing the error-weighted^ : 



2. 



X 



z 



Vsyn(*<) ~ Vdbs(Xj) 



(6) 



where the xt are the spatial positions of the observed and the 
simulated RC with veloc i ties V b s (x,) and uncertainties V eT t(Xi). 
Similar to Piz agno et alJ (120071) we add 5 kms -1 in quadrature 
to all velocity uncertainties to avoid the fit being dominated by 
high S/N data points. This is motivated by non-circular motions 
of the same magnitude, present in all disk galaxies. We use 
the photometric scale length rj and a typical value of lOOkm/s 
as initial guesses for the free parameters r t and V max . Fig. [3] 
gives an example of a best-fitting synthetic RC and the inferred 
intrinsic rotation velocities. Note that while the observed data 
points and the best-fitting RC are plotted against the pixel 
position along the spatial axis (i.e. perpendicular to the direction 
of dispersion), the intrinsic RC shows the intrinsic rotational 
velocities (see Eq.[4]i as a function of galacto-centric radius. 
Rotation curves that lack data points in the flat region tend to 
result in overestimated turnover radii and V max -values. To tackle 
this issue we constrain r, in the minimisation process (a python 
implementation of a Sequential Leas t SQuares P rogramming 
(SLSQP) optimisation, originally by iKrafJ dl988h ). We allow 
the best-fit r,-values to vary only within a cr,-interval around the 
initial value rj. We determine the scatter cr, solely using rotation 
curves that result in best-fit V max -values changing by less than 
5% when fitted with or without r, as a free parameter, cr, is 
then set to be the standard deviation of the relative differences 
between rj (i.e. r, held fixed) and r, (i.e. r, free parameter), 
multiplied by r</. We obtain a value of cr, = 0.176 * rj. 
We estimate the errors of the best-fit parameters V max and r t by 
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Fig. 3. Example of rotation-curve fitting: The blue squares indi- 
cate the observed rotation velocities as a function of the posi- 
tion along the spatial axis of the spectrum (in pix). The black 
line depicts the best-fit synthetic RC to these data points. From 
the best-fit parameters (here V max = 118km/s, r, — 4.1 pix) we 
can infer the intrinsic shape of the rotation curve (in red). Here 
the x-axis is the galactocentric radius (in kpc) (here 6 = -36°5, 
i = 68?5, z = 0.17, FWHM = 077). Geometric and observational 
effects lower the observed rotation velocities. 




0.4 0.6 

rms, 



Fig. 4. The histograms show the rms„ distribution for the low- 
(thick line) and high-quality (thin line, hatched) TF objects of 
our whole sample. For each histogram bin, the relative fraction 
of high-quality objects is depicted (black squares). It is well cor- 
related with the rms„ value of the RC fit. 



generating bootstrap samples of a given observed rotation curve 
and calculating the standard deviation for the distribution of the 
best-fit parameters obtained from each iteration. To this end, 
for each data point (at a certain position x,) the velocity value 
V bs(xi) is drawn from a Gaussian distribution with mean value 
V bs(xi) and dispersion V e „(xj). Additionally, in each bootstrap 
iteration two data points are randomly excluded from the fitting 
process to account for a possible strong influence of single 
values on the best-fit parameters. Furthermore, the inclination 
and position angle are allowed to vary for each bootstrap 
iteration according to a normal distribution around the preset 
value (with a scatter cr = 2°). This accounts for the uncertainty 
in the derivation of the inclination angle and a possible slight 
misalignment between slit direction and apparent major axis. 



3.5. Rotation-curve quality assessment 

Prior to any RC fitting, we inspect each RC visually and classify 
it according to the following scheme: non-rotating, low-quality 
Tully-Fisher object, high-quality Tully-Fisher object. The non- 
rotating class holds galaxies for which a clear distinction be- 
tween an approaching and receding side is not feasible. Only 
five galaxies fall into this class (all of them cluster members). We 
exclude these objects from further analysis. Tully-Fisher objects 
(176 in total) show reasonable disk rotation. We assign galaxies 
to the low-quality class when they have a high degree of distor- 
tion, an RC with small extent or rigid-body rotation. 
A more quantitative measure for the suitability of a galaxy as a 
TF object is provided by the (error-weighted) root mean square 
(rms) of the RC fit, i.e. the normalised sum of the squared differ- 
ences between the data points and the best-fit curve. Rotation 
curves deviating from a Courteau model (see Eq. 3} will in- 
evitably produce higher rms values than undistorted ones. Since 
absolute rms values by construction are higher for fast rotators, 



a normalised rms,, = rms/V max is a more reliable measure. Of 
course, symmetric RCs that do not sufficiently trace the flat re- 
gion or are identified as rigid body rotators, might also result in 
low rms values, despite being classified as low-quality objects. 
Hence, the correlation between the visual classification and the 
rms,, measure is not expected to be one-to-one. Fig. |4]plots the 
fraction of high-quality objects as a function of the rms,, values. 
Both measures are well correlated and result in a highly signifi- 
cant (p < 0.01%) Spearman rank-correlation of p - -0.52. 
Table Q] gives an overview of the mean rms,, values for various 
subsamples and confirms that high-quality objects have lower 
rms,, values than low-quality objects. This indicates that the lat- 
ter on average have more distorted kinematics, i.e. their devia- 
tion from an ideal Courteau model is stronger. This higher degree 
of kinematic distortion is also reflected in the mean values of the 
asymmetry measures defined in Paper I. There, w e introduced 
an asymmetry index A following Dale et al. (2001), which com- 
pares corresponding rotation velocities between the approaching 
and receding side of a rotation curve, and a visual asymmetry 
parameter A v ; sua i2, which distinguishes between three classes: 
undistorted (0), slightly distorted (1) and heavily distorted (2). 
High-quality objects (A - 14 + 1%, A visua i2 = 0.51 ± 0.06) have 
significantly lower values than low-quality objects (A = 49+6%, 
^visuai2 = 1-10 + 0.10) in both classification schemes. 

3.6. Linear regression 

Various methods are used in the literature for fitting a regression 
line to data points. We obtained a TF relation of the form: 



M - a * log V + b 



(7) 



by following iPizagno et alJ (120071) and adopting a maximum 
likelihood (ML) method to dete r mine the slop e a and intercept 
b (see also iBamford et al.ll200"6t IWeiner etakl 12006). Applying 
Bayes' theorem allows to calculate the probability of a model 
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Population 


Member 
galaxies N 


high-quality 
fraction 


rms„ 
[%] 


Total 


181 


0.62 


20.0 ±1.7 


Cluster (all) 


95 


0.58 


22.2 ± 2.5 


Cluster (hqTF) 


55 


1 


11.5 ±0.8 


Cluster (lqTF) 


35 





34.7 ±4.9 


Field (all) 


86 


0.66 


19.8 ±2.5 


Field (hqTF) 


57 


1 


11.4 ±1.1 


Field (lqTF) 


29 





36.4 ±6.3 



Table 1. Mean rms„ values from RC fitting for our sample, as 
well as for certain subsets. Cluster and field galaxies are further 
subdivided into a low and a high-quality TF class. Five cluster 
objects are classified as non-rotating and not used in the TF anal- 
ysis. 



given the observed data: 



P(model | data) oc P(data | model) x P(model), 



(8) 



where P(model) i s the prior. Flat priors are sufficient in most 
cases (D' Agostini 2005), hence they drop out of the equation. In 
contrast to ordinary least square (OLS) methods this Bayesian 
approach allows to include observational errors on both axes 
as well as an intrinsic scatter <x int . Even if measurement errors 
were negligible, a relationship between two physical variables 
will have some scatter. This intrinsic scatter reflects variations 
in physical parameter s that are no t explained by the scatter of 
the observables alone dKellvll20C)7h . We assume Gaussian distri- 
butions for the (uncorrected) observational errors as well as for 
the intrinsic scatter. We further assume a uniform intrinsic distri- 
bution of the independ ent variable x,, so th at the log-likelihood 
to maximise reads (e.g. Pizag no et al.ll2007l) : 



In L v (a, b, a- mt ) = - ^ In (of nt + a 2 cr 2 xi + cr$J) 



■I 



(y,- - axj - b) 
o 2 . + a 2 cr 2 ■ + cr 2 

int x,i y, 



(9) 



+ const.. 



where <r v , and cr xi are the measurement errors of the dependent 
variable y,- and independent variable x,, respectively. Compared 
to traditional least square methods, this fitting approach has the 
advantage that data points with large errors in either of the vari- 
ables will have less influence on the final result. Moreover, the 
inclusion of the intrinsic scatter prevents data points with very 
small obser vational er rors from dominating the fit. 
However, Kelly J2007I) emphasised that adopting a uniform dis- 
tribution for the independent variable can produce a bias in the 
ML estimate and is not appropriate for most situations (as is e.g. 
the case for the distribution of the absolute B-band magnitudes 
of our sample, see Fig. |5). Instead, they proposed to approxi- 
mate the intrinsic distribution of the independent variable using 
a mixture of Gaussian functions. For simplicity, we adopt a sin- 
gle Gaussian with mean /i and variance t 2 determined from a fit 
to the histogram of the independent variable of our whole sam- 



ple. Then, the log-likelihood of Eq.|9]becomes: 

(yi-E) 2 



In 



L G (a, b, cr mt ) = - J] In (0 2 ) - £ W Q2 + const., 



xfiT 2 Hacr 2 . 



E = b + 



©2 =fl 2 T 2 + (r 2 2 



(10) 



aV 



x,i 

where E and 2 are the expectation value and variance of y,- at 
Xt, drawn from a Gaussian distribution with mean fi and variance 
t 2 . The covariance between the observational errors a- y j and cr x j 
was set to zero. For t 2 — > oo, i.e. an infinitely broad Gaussian, 
Lq of Eq. PTOlconverges to Ljj of Eq.|9] 

We are aware that the intrinsic distribution of luminosities would 
rather be a Schechter-like function. However, a Gaussian is a 
good approximatio n to the observed lu minosity function of spi- 
ral galaxies only (Bingg eli et al.| [l988). Furthermore, it is in any 
case more realistic than a flat distribution. 
We perform the minimisation of the (negative) log-likelihood us- 
ing a downhill simplex algorithm. We estimate the errors of the 
best-fit parameters a, b and <x int utilising bootstrap simulations 
(« = 1000). 

Instead of fitting the "forward" TFR of Eq. [7] i.e. y = M and 
x - log V, many studies use the "inverse" relation: 



1 b 

log V = -M - - 

a a 



(11) 



The inverse fit, i.e. the velocity treated as the dependent vari- 
able, is less prone to incomple teness biases that arise f rom the 
magnitude limit of a sample dTullv & Courtoisl 1201 2l IWillick 
1994). Kellv (2007) demonstrated that if the sample selection 
is based on the independent variable of the fit (the luminosity 
in our case), then the best-fit parameters are not influenced by 
selection effects. However, treating either the luminosity (for- 
ward fit) or the velocity (inverse fit) as the dependent variable 
does not reflect the situation properly, since by construction in 
each case the scatter in the dependent variable is minimised. 
To study the underlying relationship between two observables 
a method treating them symmetrically could be more appro- 
priate. Such methods utilise both the forward and the inverse 
best-fit parameters to calculate a best-fit slope that lies between 
the two extremes. Ilsobe et all dl99 0) describe for instance bi- 
sector regr ession, orthogona l regre ssion or reduced major-axis 
regression. iTullv & Courtois (2012) simulated the selection bias 
for different fitting methods and recommend the inverse relation 
when using the TFR as a means of distance measurement. On the 
other hand, to study the physical impacts on the TFR a bivariate 
fit would be a good compromise. 

We choose to use the inverse fit when comparing different popu- 
lations or distant to local studies. However, for the investigation 
of Tully-Fisher residuals (luminosity offsets) we use the bisec- 
tor parameters cal culated from the forward and inverse fits as in 
Ilsobe et all d!990l) . Whatever the type of fit, we always report the 
resulting parameters in terms of a forward relation. 

3. 7. Incompleteness bias correction 

By construction magnitude-limited samples result in an incom- 
plete coverage of the faint end of the luminosity function. In a 
TF sample with magnitude limit faint slow rotators are under- 
represented. In the end, this leads to an underestimation of the 
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Fig. 5. B-band luminosity histogram of our spectroscopic sam- 
ple. The red solid line depicts a Hanning-smoothed fit to this 
histogram. The blue dashed line indicates a Schechter fit to the 
bright end [Mb < — 19?5). The obtained parameter is M* - 
-20^9 for the characteristic luminosity. We assume a slope at 
the faint end of a — -1.00. The incompleteness of our sample is 
given by the difference between the two lines. 



TFR slope accompanie d by an overestimation of the intercept. 
iGiovanelli et al.l (1 19971) developed a method to account for this 
bias. First of all, we have to assume an intrinsic luminosity dis- 
tribution according to a Schechter luminosity function. Fig. [5] 
shows a histogram of the B-band luminosities and a Hanning- 
smoothed fit to this histogram. The latter is then compared to a 
fit of a Schechter luminosity function to the bright end of the data 
(Mb < — 19^5). For these luminosities, the sample is assumed to 
be complete. The ratio between the two lines (maximum set to 
1) provides us with a measure c(Mb) for the completeness of the 
sample. 

To simulate the incompleteness bias we generate Monte-Carlo 
simulations with 10 4 iterations. For a given rotation velocity a 
luminosity is calculated according to an intrinsic TFR. The TFR 
scatter is simulated by randomly drawing from a normal distribu- 
tion with a dispersion equal to the scatter in our data. The lumi- 
nosity value Mb thus obtained now has the probability c(Mb) to 
enter the sample. Of course, fainter galaxies have a lower chance 
to be considered. The average residual from the TFR gives an es- 
timate of the magnitude bias for a given rotational velocity (see 
Fig.|6]l. The magnitude bias AMb - /(V2.2) is then approximated 
by a linear fit. To account for the slightly different luminosity and 
redshift distributions of the field and cluster sample, the magni- 
tude bias was computed for the two populations separately. 



4. Tully-Fisher relations 

Quantitative comparisons between different Tully-Fisher studies 
are very difficult. First of all, there can be differences in the ob- 
servational approach. The rotation of a galaxy may be traced e.g. 
by radio data of the H I disk or by optical emission lines stem- 
ming from HII regions. Radio observations have the advantage 
of tracing the rotation of a galaxy to large radii, where the dar k 
matter halo dominates the mass density (e.g. lVerheiienl l2001). 
Pizagno et al. (2007) pointed out that obtaining velocity-widths 
from H I rotation curves results in rotation velocities systemat- 
ically higher for fainter galaxies compared to optical observa- 
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Fig. 6. Simulated incompleteness bias obtained from Monte- 
Carlo simulations. The data points depict the average TF residual 
for a certain velocity. The blue solid line indicates a linear fit to 
these data points. 



tions. This leads to a shallower TFR. Furthermore, some stud- 
ies use integrated line-width while others exploit spatially re- 
solved rotation curves or even velocity fields. Compared to the 
traditional long-slit spectroscopy, IFU (Integral Field Unit) ve- 
locity fields might be superior in identifying galaxies with dis- 
torted and peculiar kinematics. Besides this, differences in sam- 
ple selection (magnitude limit, morphological type, constraints 
on kinematic distortion), different fitting methods and consider- 
ations of observational effects are the most common sources of 
discrepancies. To avoid a possible bias, we focus on the compar- 
ison of subsamples in our analysis, e.g. between field and cluster 
galaxies or the blue cloud and dusty red population. 
Studies designed to use the TFR as a distance indicator (usually 
with the goal of constraining the Hubble constant H$), need a 
pruned sample containing only "model" galaxies such as undis- 
torted late-type spirals. This guarantees a tight TFR with a small 
scatter. However, we also intend to investigate galaxies that de- 
viate from this ideal, which might reveal information on interac- 
tion mechanisms in the cluster environment. 

4. 1 . Basic fits 

Table [2] gives an overview of the parameter estimates using the 
maximum likelihood method with a Gaussian intrinsic distribu- 
tion of the independent variable (Eq.fTOt. Note that using a uni- 
form intrinsic distribution instead (Eq. not shown in the table) 
yields consistent values within the lcr errors. 
First, we concentrate on the B-band TFR using only rotation 
curves of our high quality class. Inverse fits to the sample of 
field (57 objects, a = -7.51 + 0.81) and cluster (55 objects, 
a - -6.58 + 0.61) galaxies yield consistent slopes within the 
lcr-errors. Hence, in the following subsections we mainly use a 
fixed slope, determined from an inverse fit to the combined sam- 
ple of cluster and field galaxies. This is a common approach to 
reduce the number of free parameters and is justified by the fact 
that most other TF studies also do not find signifi cant differences 
in the slope between cluster and field galaxies (Bamford et al 
2006; Moranetal. 2007; Nakamura et al. 2006; Ziegleretal. 



2003). Fig.|7]shows an inverse fit to the high-quality objects of 
our whole sample. As a comparison we also depict the (shal- 
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B-band TFR 



Stellar-mass TFR 



slope a 



intercept b scatter (cr int ) 



slope a intercept b scatter (<x int ) 



Total (hq, 112 galaxies) 

forward -5.10 ±0.30 -8.87 ±0.61 0.66 ±0.05 

inverse -7.04 ± 0.46 -4.95 ±0.95 0.77 ±0.10 

bisector -5.92 ±0.38 -7.20 ±0.77 0.71 ±0.07 
Cluster (hq, 55 galaxies) 

forward -5.28 ±0.32 -8.33 ±0.64 0.57 ±0.07 

inverse -6.58 ± 0.61 -5.69 ±1.25 0.64 ±0.10 

bisector -5.86 ±0.45 -7.12 ±0.93 0.60 ±0.08 
Field (hq, 57 galaxies) 

forward -5.03 ± 0.44 -9.15 ± 0.89 0.69 ± 0.07 

inverse -7.51 ±0.81 -4.20 ±1.66 0.84 ±0.15 

bisector -6.03 ± 0.61 -7.17 ± 1.22 0.76 ±0.11 
Total (h+lq, 176 galaxies) 

forward -4.09 ± 0.38 -10.85 ±0.74 0.90 ±0.08 

inverse -7.08 ± 0.64 -5.03 ± 1.30 1.20 ±0.17 

bisector -5.20 ±0.50 -8.96 ±0.97 1.04 ±0.12 



2.41 ± 0. 12 4.95 ± 0.24 0.26 ± 0.02 

3.04 ±0.19 3.69 ±0.39 0.29 ±0.03 

2.69 ±0.15 4.40 ±0.31 0.27 ±0.03 

2.34 ±0.16 5.14 ±0.33 0.25 ±0.02 

2.86 ± 0.21 4.07 ± 0.42 0.27 ± 0.03 

2.57 ±0.18 4.65 ±0.38 0.26 ±0.03 

2.48 ±0.16 4.81 ±0.34 0.26 ±0.04 
3.25 ±0.36 3.27 ±0.75 0.30 ±0.06 
2.82 ±0.25 4.14 ±0.50 0.28 ±0.05 

1.94 ±0.16 5.95 ±0.32 0.43 ±0.04 

3.37 ±0.35 3.15 ±0.70 0.58 ±0.09 

2.49 ±0.25 5.01 ±0.48 0.50 ±0.06 



Table 2. Parameter estimates using the maximum likelihood method with a Gaussian intrinsic distribution of the independent vari- 
able. For each subpopulation forward, inverse and bisector parameters are listed. All values are given in the forward representation 
of the Tully-Fisher relation. 



lower) slopes of the bisector and forward fit as well as the posi- 
tion of the low-quality galaxies. 

Note that an inverse fit minimises the scatter with respect to the 
rotation velocity and not to the luminosity. Hence, an intrinsic 
scatter obtained from an inverse fit and transformed into a value 
with respect to the luminosity is always larger than if directly 
obtained using a forward fit. This needs to be considered when 
comparing to studies that adopt a forward fit. Furthermore, re- 
ducing the number of free parameters (e.g. holding the slope 
fixed) can also slightly change the best-fit intrinsic scatter. 
For the stellar-mass TFR the parameters of field and cluster 
galaxies also agree within the error bars. 

Including also the low-quality objects in the fit (176 objects in 
total) increases the intrinsic scatter significantly (see Table |2j. 
Interestingly, the best-fit (inverse) values for slope and intercept 
barely change. This is mainly due to the weighting capability of 
our linear regression method. The adopted bootstrapping method 
for RC fitting (see end of Section 13.4b assigns larger errors to 
low-quality objects {6V = 19 ±3 km/s) compared to high-quality 
objects (6V - 10 ± 1 km/s). 

For fixed slopes (a = -7.04 and a - 3.04 for the B-band and 
stellar-mass TFR) we find a difference of AM B = -0 m 70 ± 0'. n 28 
and AM, = 0.46 ± 0.16 when comparing the TF intercept be- 
tween low- and high-quality objects individually, i.e. distorted 
objects are preferentially located above the TFR of high-quality 
objects. This means they are on average too luminous/stellar- 
massive for a given velocity or have too low a velocity for a 
given luminosity/stellar mass. The comparable shifts (AV2.2 ~ 
0.12±0.5, when transformed into a velocity shift) for the stellar- 
mass and the B-band TFR as well as a similar increase in the 
intrinsic scatter (cr int ~ 0.18 ± 0.3 in velocity) indicate that 
mainly an underestimated rotation velocity causes this shift for 
low-quality objects. 



4.2. Cluster vs. field relation 

Holding the TFR slope fixed at a — -7.04, we now compare the 
B-band TFR of high-quality cluster and field galaxies. We find 
an offset between the two populations (Fig. [8] left panel) with 
field galaxies being slightly brighter than their cluster counter- 
parts (AMfl = -0^42 ± fT15). 

Furthermore, we find a slight difference between the mean SED 
types. Field galaxies have on average a bluer SED than cluster 
galaxies ((r type )=5.8 vs. 5.2 ±0.2), which could explain part of 
the offset in the TFR. T-types of 5 and 6 correspond to typical 
spectra of Sc and Scd spirals, respectively. A linear fit to the rela- 
tion of the B-band residuals as a function of the SED-type yields 
a brightening of AM B ~ -OPQ9 ± 0^05 for an SED-type change 

of Ar tyP e = 1. 

Another aspect might be the deviating redshift distribution 
of both samples. According to a trend seen in previous TF 
evolution studies (e.g. Bamford et al. 2006; Bohmetal. 200' 
iFernandez Lorenzo et al.: 2010: Mi ller et alJl201 it IWeiner et 



2006; Ziegler etal] l 2002T 



r- 

2 



the slightly higher redshift of our 
field galaxies ((z)=0.245 vs. 0.165) coincides with higher star- 
formation rates (e.g. Hopkins 2004) and consequently higher B- 
band luminosities. Most studies measure a brightening in the 
B-band consistent with AM B - -1™0 ± 0^5 from z - to 
z ~ 1. Galaxy form ation models based on ACDM cosmology 
(lDuttonetalJl201lh produced similar values {AM B - -07 1 9). 
Hence we would expect an offset of AM B ~ -O^l ± 07*1 stem- 
ming from the redshift difference between our cluster and field 
sample. 

Altogether, systematic differences in the composition of the clus- 
ter and field sample may account for half of the observed offset 
in the B-band intercept. The cluster environment associated with 
a quenching of star-formation could be responsible for the re- 
maining difference (AM B ~ -0^23 ± (T19). Of course this dif- 
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Fig. 8. TFR for the sample of high-quality cluster (55 objects; blue squares) and field (57 objects; green triangles) galaxies. Open 
symbols indicate the corresponding low quality objects. The blue dashed and green dot-dashed lines indicate the best-fit inverse 
regression lines of cluster and field galaxies, respectively. B-band TFR (left panel): Field galaxies are brighter by AMb - d"42 + 
0^15. SmTFR (right panel): We find no difference between the intercept and intrinsic scatter of the two populations (AM, = 
0.00 + 0.07, Acr int = 0.01 + 0.04). 



ference is barely significant. Furthermore, if the assumption of 
a fixed slope between cluster and field galaxies is not fulfilled, 
i.e. cluster galaxies have an intrinsically shallower slope, then an 
underestimation of the TF intercept is expected, when fixing the 
slope. 

On the other hand, under the assumption of such a luminosity- 
redshift correlation, the broader redshift distribution of the field 
sample (<x z =0.048 vs. 0.005 for the cluster, which is a factor 
of ~10) should result in a slightly larger TF scatter. The intrin- 
sic scatter amounts to <x int = 0.68 ± 0.06 for the cluster and 
o"int = 0.79 + 0.08 for the field population, which is in agreement 
within the errors. Given the redshift distribution of our field sam- 
ple and assuming an evolution of AMb = (-1™0 + 0. n 5) x z (e.g. 
Bamford et al. 2006), simulations yield an estimate of Act = 
0.05 + 0.02 for the increase of the scatter due to the broader 
redshift distribution. 

However, the situation changes when only low-quality Tully- 
Fisher objects are considered. While the difference in the in- 
tercept is similar (AM B - -0. n 62 + 0? 1 48), the intrinsic scatter 
for the cluster sample (35 objects; o- m - 2.20 + 0.32) is sig- 
nificantly higher than for the corresponding field subsample (29 
objects; <x int = 1 .42 + 0.25). Whereas cluster galaxies with more 
regular kinematics show no enhanced scatter, cluster-specific in- 
teraction processes seem to induce a larger TF scatter in galaxies 
with distorted kinematics. 

However, when we examine the high-quality stellar-mass TFR, 
the differences between cluster and field galaxies disappear. 
Again, we fix the slope to the (inverse) value of the combined 
sample (a = 3.04). Then, the intercept (AM, = 0.00 ± 0.07) and 
intrinsic scatter (Acr int = 0.01 + 0.04) for both populations are 
in total agreement. Fig. [8] (right panel) illustrates the coinciding 
sample of cluster and field galaxies for the stellar-mass TFR. 
This finding suggests that any measured differences between the 
TFR of cluster and field galaxies stem from different M/L ra- 



tios and hence star-formation rates affecting the B-band lumi- 
nosity. It confirms that differences in luminosity and not in ve- 
locity induce the offset of the B-band TFR between cluster and 
field galaxies. 



4.3. Local comparison samples 

In the following we want to compare our sample to local TF stud- 
ies. We select th e B-band data of Tully & Pierce (2000, Sample 
T hereafter) and Verheijen (2001, Samp l e V h ereafter) as well 
as the stellar-mass data of Pizagn oet al.l (120051 Sample P here- 
after) as comparison samples. Both B-band studies exploited 
H I linewidth measurements W20 (at 20% of the peak flux) to 
derive rotation velocities, and applied the same intrinsic mass- 
depen dent absorption correction as we do. The Tull y & Piercel 
(2000) sample comprises 115 spirals (Sa-Sd) in various envi- 
ronments and was mainly designed to calibrate the TFR for 
subsequent Ho determination. We also use the total sample of 
45 galaxies from lVerheijenl (120011) as a less pruned comparison 
sample, including ob jects with potentially distorted kinematics. 
Piza gno et al.l (2005) selected 81 disk-dominated galaxies from 
the SDSS and derived the rotation velocity at a radius of 2.2 disk 
scale lengths, V2.2, from H a rotation curves. 
By applying the same fitting procedure to our sample and the 
local samples, we rule out a potential source of discrepancies. 
Table[3]lists the TF parameters we obtain when performing a free 
inverse fit and when fixing the slope to a=-7.04 (B-band TFR) 
or a=3.04 (smTFR): Fig.|9]shows an overview of the data. Both 
local B-band samples (left panel) have a fainter intercept than 
our sample. Fixing the slope yields differences in the intercept 
of AM B = 0^33 + (T 10 (Sample T) and AM B = 0M3 + 0?13 
(Sample V), respectively. This, as well as the smaller scatter, 
can be most likel y attributed to the hig her redshift of our sam- 
ple. Compared to iTullv & Piercel (l2000t) the less strict selection 
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Fig. 9. TF data for our total sample (blue points, solid line) and the local comparison samples. The depicted regression lines are for 
the fixed slope of our sample. B-band (left panel): the intercept differences ar e AM B - (F33 ± 07*10 for the JTullv & Piercel (feOOO, 
black triangles, dashed line) sample and AMb - 0P43 + 0^13 for the sample of Verheiien (2001 , red circles, dot-dashed line). Stellar 
mass (right panel): the intercept of our sample agrees with the local sample of Pizag no et alj (12005, black triangles, dashed line) 
(AM, = 0.01 ± 0.04). 



Local TFR 



slope a intercept b scatter (cr int ) 



Sample T 



Sample V 



Sample P 



B-band TFR 

-7.93 + 0.29 

-7.04 fixed 

B-band TFR 

-6.98 + 0.52 

-7.04 fixed 

stellar-mass TFR 

3.40 + 0.11 

3.04 fixed 



-2.66 + 0.64 
-4.62 + 0.04 

-4.65 + 1.11 
-4.52 ± 0.09 

2.89 ± 0.23 
3.70 ± 0.02 



0.43 ± 0.04 
0.40 + 0.03 

0.61 ±0.10 
0.61 ± 0.09 

0.17 ±0.02 
0.15 ±0.02 



Table 3. Parameter estimates of local comparison samples using 
the maximum likelihood method with a Gaussian intrinsic dis- 
tribution of the independent variable. For each sample we carry 
out an inverse fit with a free and a fixed slope. 



constraints might also contribute to the larger scatter in our data. 
The best-fit slopes are in good agreement with our values, al- 
though the relation of the Sample T is slightly steeper. A lack of 
faint fast-rotating galaxies in our sample (especially in the clus- 
ter population) seems to be responsible for this difference. As 
already mentioned this could (in contrast to velocity width com- 
puted from H I data) partly be due to t he derivation of rota tional 
velocities from optical emission lines (Pizag no et al.l l2007). 
The right panel of Fig. [9] sho ws the stel l ar-m ass TFR in compar- 
ison to the local data of Pizagno et al. (2005). When fixing the 
slope, the difference in the intercept between the local P sam- 
ple and our sample is negligible (AM* = 0.01 ± 0.04). This is 
in agreement with previous studies that find no significant evo- 



lution of the stellar-ma ss TFR with redshift up to z < 2 (e.g. 
lMilleretal.ll20TTll2012h . 



4.4. Blue cloud vs. dusty red galaxies 

In the U-V colour-magnitude diagram, dusty red galaxies over- 
lap with the red sequence of ellipticals and lenticulars without 
star formation. This is due to significant dust extinction, similar 
to blue cloud galaxies, along with a low-level of star formation 
(Wolf et al. 2005). They are mainly early-type spirals and their 
star-formation is ~four times lower than in blue cloud spirals of 
similar mass (I Wolf et alJl2009t) . Hence, we expect that dusty red 
galaxies are located below the TFR of blue cloud galaxies. Fig. 
[10] (left panel) confirms this expectation. The dusty red galaxies 
preferentially populate the region below the (bisector) regres- 
sion line of the whole high-quality cluster sample (blue clouds 
+ dusty reds). The mean (error-weighted) residual amounts to 
AM B = 0? 1 24 ±0^11 compared to AM B = -0?11 ± <J?06 
for blue cloud galaxies, i.e. dusty red galaxies are on average 
AMb = 0^35 ± (J?\2 fainter for a given velocity than blue cloud 
galaxies. The trend that the cluster slope is slightly shallower 
might stem from the high fraction of dusty red galaxies in our 
sample. They populate the intermediate-to-high mass regime of 
our sample and lie on average below the TFR (see FigJTO] left 
panel), thus reducing the slope. Excluding the dusty red galaxies 
(12) from the fit yields a slope closer to the field value (43 ob- 
jects, e.g. for the inverse fit a — -7.07 ± 0.74). 
The right panel of Fig. [TOl shows the same comparison between 
dusty red and blue cloud galaxies, but this time for the stellar- 
mass TFR. Here, the mean (error-weighted) residuals from the 
bisector regression line are AM* = -0.05 ± 0.05 for dusty red 
and AM* = 0.01 ± 0.02 for blue cloud galaxies, respectively. In 
contrast to the B-band residuals, this is an insignificant differ- 
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5. Tully-Fisher residuals 
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Fig. 7. B-band TFR for the combined sample of cluster (55) and 
field (57) galaxies. Only high-quality objects are considered. The 
blue solid line indicates the best-fit inverse regression line. The 
blue dotted lines depict the corresponding \cr and 3<x intervals, 
respectively. The green and red line are the best-fit forward and 
bisector fits, respectively. Black open circles indicate the posi- 
tions of the low-quality objects. 



ence (AM* = 0.06 + 0.05). This suggests that at given rotational 
velocity dusty red galaxies have similar stellar masses as blue 
cloud galaxies. A difference only appears in the B-band TFR as 
a result of different M/L values. 

The offset in the average B-band residuals fits into the scenario 
that dusty red galaxies might be the progenitors of SO galax- 
ies. The latter are found to systema tically lie below the Tully - 
Fisher relation of spiral galaxies (Aragon-Salamanca 2008). 
Cortesi et al. (2013) investigated the kinematics of six SO galax- 
ies exploiting spectra of planetary nebulae. Besides a position of 
the disk components below the TF relation of local spirals, they 
also found that the bulge components populate the region above 
the Faber-Jackson relation of elliptical galaxies. In other words, 
bulges are too bright but the disks are too faint for a given ve- 
locity measure. Furthermore, the disks exhibit an increased frac- 
tion of random motion compared to spiral galaxies. The authors 
pointed out that neither a quiescent scenario like disk-fading nor 
a violent interaction like merging are likely explanations of the 
observations. A milder process like ram-pressure would, how- 
ever, be well in compliance with what we find. 
Our sample of 182 disk galaxies contains 14 objects (11 in the 
cluster; 3 in the field) with an SO morphology. 22% of the cluster 
galaxies with spiral morphology are dusty red galaxies, and 9 of 
1 1 lenticulars (82%) have a dusty red SED-type. Only four clus- 
ter SO galaxies show high-quality rotation curves. Their mean 
(error-weighted) residual is AM# = 07*74 ± 0^23 which is even 
fainter than the average dusty -red population. The small number 
statistics do not allow to draw conclusions here, but the numbers 
fit into the transformation scenario. 



Several effects might influence the position of a galaxy in the TF 
plane. On the one hand, enhanced star- formation or higher dust- 
extinction rates can shift a galaxy along the luminosity axis, but 
do not affect stellar mass (changes in M/L and L cancel out). 
On the other hand, a high fraction of non-circular motion, pos- 
sibly caused by kinematic distortions, can change the radial ve- 
locity profile and the measured rotational velocity. Variations in 
the stellar mass-to-light ratio or stellar-to-total mass ratio induce 
shifts in the TF plane (Kannappan et al. 2002). 
In this section we mainly focus on the analysis of the B-band 
AMb and stellar-mass residuals AM„. We examine their correla- 
tion with other parameters and a possible dependence on envi- 
ronment. We focus on the population of dusty red galaxies and 
investigate whether differences to blue cloud galaxies exist. As a 
proxy for environment we adopt the scaled cluster-centric radius 
r s - '"/'"200 from Paper I. For a given subcluster, r is the projected 
cluster-centric distance and r2oo is the radius where the mean in- 
terior density is 200 times the critical density. 
Furthermore, we also wish to have a measure for the asymme- 
try of the stellar disk at hand and we adopt the index A morp h, 
which is computed by comp aring an image to itself rotated 
by 180° (Cons elice et al.ll2000b . A morp h quantifies the degree of 
morphological distortions (see also Paper I). Fig. [TT] shows the 
normalised rms as a function of the morphological asymmetry 
A m orph- This plot resembles a plot from Paper I (Fig. 13 + 17), 
but this time the rotational asymmetry is parametrised by rms„, 
i.e. the deviation of an observed rotation curve from its best-fit 
Courteau model. More kinematically distorted field galaxies ex- 
hibit higher morphological asymmetry. A Spearman rank-order 
correlation test yields a significant (p - 2.7%) correlation of 
p = 0.31. On the other hand, dusty red galaxies (p = -0.58, 
p - 4.4%) have higher rms„ for a low A morp h. Although blue 
cloud galaxies show no significant correlation, they have - com- 
pared to field galaxies - also increased rms„ values for smooth 
stellar disks. Following our argument in Paper I, interaction pro- 
cesses that affect the morphology (galaxy-galaxy interactions, 
harassment) most likely also distort the gaseous disk of a galaxy 
and hence perturb its rotation curve. Depending on the ICM 
density and the rela tive velocities, ram-pressure (oc pic M * v*, 
iGunn & Gottl 119721) is able to significantly influence the gas 
content of a galaxy and strip part of its material, while not af- 
fecting its stellar disk (IChunget all 120 08: Ken nev et all l2004t 
iKronberger eTall 12008b; Ouilis et al. 200CJ). In the field, the 
rms n of the RC-fitting positively correlates with morphologi- 
cal disturbances. There, low relative velocities of galaxies favour 
tidal interactions. For morphologically undistorted cluster galax- 
ies and especially dusty red galaxies a cluster-specific mecha- 
nism that effectively distorts the rotation curve but leaves the 
stellar disk unaffected is required. Ram-pressure stripping is the 
best candidate process. 

Interestingly, low-quality and high-quality objects have the 
same mean morphological asymmetry (Amcphjq = 0. 20 + 0.01; 
Amorph .ha = 0.22 + 0.02). This is in agreement with Jaff e et all 
(201 1), who also do not find a direct link between the morpho- 
logical asymmetries in the stellar disk and the fraction of good 
rotation-curve fits. 

This raises the question of whether the stellar asymmetry of a 
galaxy influences its offset from the best-fit TFR at all. Fig. 
[12] (left panel) shows that for the class of high-quality objects 
a higher morphological asymmetry induces a brighter B-band 
residual. The correlations are highly significant (p < 0.1%) for 
the field (p = -0.47) and the cluster blue cloud (p - -0.45) 
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Fig. 10. TFR for the sample of high-quality cluster galaxies, separated into dusty red (red stars) and blue cloud (blue squares) 
galaxies. The black line indicates the best-fit bisector regression line for the whole cluster sample. Fixing the slope to this value, 
the red dashed-dotted and blue dashed lines depict the corresponding best-fit lines for dusty reds and blue clouds, respectively. 
B-band (left panel): For a given velocity dusty red galaxies are on average AMg = CP35 ± (1^12 fainter than blue cloud galaxies. 
Stellar mass (right panel): The difference between the average residuals of dusty red and blue cloud galaxies is not significant 
(AM, = 0.06 + 0.05). 
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Fig. 12. Mean residuals of the B-band TFR for the sample of high-quality cluster galaxies, separated into dusty red (red stars) and 
blue cloud galaxies (blue squares) as well as for the field sample (green triangles). Left panel: The correlation with morphological 
asymmetry is highly significant for blue cloud and field galaxies. Mid panel: Galaxies with high specific star-formation rates tend 
to lie above the TFR. Right panel: Dusty red (p - -0.50, p - 4.0%) and blue cloud (p - 0.39, p - 0.5%) galaxies show opposite 
correlations as a function of cluster-centric distance. 



population. Dusty red galaxies, on the other hand, produce no 
significant correlation (p - 0.20, p - 27.9%). The same be- 
haviour is seen for low-quality objects (not shown in the figure). 
Galaxies with a higher morphological asymmetry index tend 
to have a more late-type SED and higher star- formation rates, 



which, in turn, leads to brighter residuals. Also iGiovanelli et all 
(e.g. 1997) found TF offsets slightly fainter for more early-type 
galaxies. This could explain part of these correlations. On the 
other hand, galaxy-galaxy interactions, accompanied by mor- 
phological distortions, can induce a star-formation enhancement 
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Fig. 11. The rms„ from RC fitting as a function of morphological 
asymmetry A moip h for blue cloud (43), dusty red (12) and field 
galaxies (57). Only high-quality objects are considered. 



and consequently result in larger TF offsets (e.g. lAlonso et alj 
2004; Lamba s et al]l2003l) . 

Fig. [12] (mid panel) confirms that galaxies with higher specific 
star-formation rates preferentially populate the region above 
the B-band TFR. Dusty red galaxies have reduced specific star 
formation rates. Note that we only have robust star- formation 
rates (determined from UV and infra-red data) of 30/55 high- 
quality cluster (22/43 blue clouds; 8/12 dusty reds) and 22/57 
high-quality field galax ies. Our finding is in agreement with 
iKannappan et a l. (2002) who also detect a correlation between 
H a equivalent widths (a proxy for specific star formation) and 
TF residuals. 

We discuss possible environmental effects on the B-band TF 
residuals in Fig. [12] (right panel). The mean B-band residuals 
of dusty red (p = -0.50, p = 4.0%) and blue cloud (p = 0.39, 
p - 0.5%) galaxies show a different behaviour as a function of 
cluster-centric radius. While dusty red galaxies become fainter 
towards the cluster center, blue cloud galaxies tend to be more 
luminous at small r s . An enhanced ICM density and ram pres- 
sure may trigger star formation in the gas-rich blue clo ud popula- 
tion (Kronberger et al. 2008a; Steinhauser et al. 2012), while the 
gas-poor dusty red galaxies might already be depleted of most of 
their gas due to ram-pressure stripping. 

The denser environment seems to affect the residuals of cluster 
galaxies. Galaxies closer to the cluster centre deviate more from 
the TFR, which is confirmed by the fact that the absolute B- 
band residuals of high-quality objects increase towards the clus- 
ter centre (p = -0.27, p = 1.9%). Basically, higher absolute 
residuals produce a larger scatter. Fig. [13] shows the intrinsic 
B-band scatter of high-quality objects for three cluster-centric 
radius bins. Since we detect the same trend for dusty red and 
blue cloud galaxies we only show the combined cluster sample. 
The TFR of galaxies at lower r s has a larger intrinsic scatter. 
This might be due to a mix of several effects. Gas-rich (blue 
cloud) galaxies have increased star-formation rates towards the 
cluster centre, while gas-poor (dusty red) galaxies are located 
below the TFR, possibly due to quenching in progress. The ob- 
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ple of high-quality cluster galaxies. The black dashed horizontal 
line depicts the intrinsic scatter of the total cluster sample. 
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Fig. 14. The rms„ from RC fitting as a function of cluster-centric 
radius r s for the sample of high-quality cluster galaxies. The cor- 
relation is significant (p - -0.30, p - 4.3%). 

served trends in the B-band residuals are not significant for the 
corresponding stellar-mass residuals. Hence, these trends in the 
B-band probably are mainly due to M/L- variations. A higher 
amount of RC distortions may also be responsible for the larger 
scatter at smaller r s ; galaxies residing closer to the cluster center 
have on average higher rms„ values (p = -0.30, p - 4.3%, see 
Fig-EJl, i.e. their rotation curves deviate more form a Courteau 
model. 



6. Summary and conclusions 

We continue the work of Paper I based on long-slit spec- 
troscopy of disk galaxies around the multiple cluster system 
Abell 901/902 taken with the VLT instrument VIMOS. We ex- 
ploit up to four emission lines per object to extract rotation 
curves of 96 cluster and 86 field galaxies. Assuming a tilted 
ring model and using input parameters such as the inclination 
angle, the location of the apparent major axis with respect to 
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the direction of dispersion, the width and tilt angle of the slit, 
and the FWHM of the seeing disk, we simulate long-slit spec- 
troscopy. Adopting a Courteau model for the radial rotation law 
of a galaxy, this simulation enables us to find the best-fit in- 
trinsic maximum rotation velocity given the observed rotation 
curve and the input parameters. We assess the quality of a ro- 
tation curve visually (distinguishing between high-quality and 
low-quality objects) and by utilising the rms of the rotation- 
curve fitting. We then use the rotation velocity to build a B-band 
Tully-Fisher relation (with absolute magnitudes derived from 
HST/ACS V-band images) and a stellar-mass Tully-Fisher re- 
lation. We fit a regression line adopting a maximum likelihood 
method based on Bayes' theorem. Eventually, we analyse the 
TF relations and examine B-band and stellar-mass residuals. We 
focus on differences between cluster and field environment, or 
between dusty red and blue cloud galaxies in the cluster. Dusty 
red galaxies are mainly found in the cluster environment and 
populate the red sequence in an U-V colour-magnitude diagram. 
This is partly due to significa nt dust extinction and partly due to 
reduced star formation rates (Wolf et al. 2005). The main results 
of this work are: 



1. The best-fit TF parameters (slope, intercept, intrinsic scat- 
ter) between cluster and field galaxies agree within the er- 
ror bars. However, if we assume a fixed slope, high-quality 
field galaxies are AM B - -0M2 ±07*15 brighter than the 
corresponding cluster population. Accounting for the higher 
redshift and on average bluer SED of the field population 
explains half of this effect. Since we find no offset for the 
stellar-mass TFR, the residual difference might stem from 
M/L-variations due to environmental effects. 

2. When considering low-quality objects only, the intrinsic 
scatter is increased for cluster galaxies compared to field 
galaxies. This indicates more frequent kinematic distortions 
in the cluster environment. 

3. Compared to local TFRs we obtain a higher intrinsic scat- 
ter and a slightly shallower slope. At fixed slope, we find a 
B-band offset of AM B < -OTO + O^IO. Since we find no dif- 
ference in the stellar-mass intercept, this is, at least in part, 
due to a luminosity evolution with look-back time (~ 2 Gyr). 

4. While there is no significant difference in the average stellar- 
mass residuals, dusty red galaxies are AM B - 0^35 + 0742 
fainter than blue cloud galaxies for a given velocity. This 
indicates different gas contents between the two SED-types. 

5. Exploiting the normalised rms of the RC fitting, i.e. the devi- 
ation of a RC from its best-fit Courteau model, field galaxies 
are kinematically more distorted towards higher morpholog - 
ical asymmetry. In contrast, cluster galaxies (and especially 
dusty red galaxies) with a smooth stellar disk show increased 
kinematic disturbances compared to field galaxies. A cluster- 
specific interaction mechanism like ram-pressure stripping, 
which does not affect the morphology but only the kinemat- 
ics, is the most likely explanation for this trend. 

6. Galaxies with higher specific star-formation rates have on 
average higher B-band residuals. 

7. Dusty red galaxies are fainter towards the cluster centre, 
whereas the B-band residuals of blue cloud galaxies become 
brighter for lower cluster-centric radii. While the SFR of 
gas-rich galaxies might be enhanced in denser environments, 
the gas-poor dusty red galaxies might be in the progress of 
quenching due to ram-pressure stripping when falling into 
the cluster. 



8. Galaxies residing close to the cluster centre show a larger 
TFR scatter. This is probably caused by increased/decreased 
star-formation rates and/or distorted kinematics. 

Altogether, the results of the TF analysis are consistent with the 
scenario presented in Paper I, which suggests that cluster galax- 
ies are affected by the dense environment they are residing in. 
Ram-pressure is an important mechanism especially for dusty 
red galaxies, which might be an intermediate stage in the trans- 
formation of infalling field spiral into cluster lenticular galaxies. 
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